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PRESSURE-SENSITIVE SYSTEM FOR GAS- TEMPERATURE CONTROL 
By Norman Matz and Richard S. Cesaro 


SUMMARY 

A thermodynamic relation ia derived and simplified for use 
as a temperature- limiting control equation involving measurement 
of gas temperature before combustion and gas pressures before and 
after combustion. For critical flow in the turbine nozzles of 
gas-turbine engines, the control equation is further simplified 
to require only measurements upstream of the burner. Hypothetical 
control systems are discussed to illustrate appiica C l Oil ua the 
control equations. 


INTRODUCTION 

The problem of gas- temperature control is becoming increasingly 
important for high-temperature application such as the gas-turbine 
.power plant for aircraft. The acceleration of conventional gas- 
turbine power plants is currently induced by increasing the fuel 
flow to the engine. Because fuel flow is increased before the air 
flow increases, such a process of acceleration is accompanied by 
an increase in combustion-gas temperature to values that may exceed 
the temperature limitation of the engine. Temperature- limiting 
control is therefore required in order to eliminate the possibility 
of engine failure from overheating during this transient phase of 
operation. 

The temperature limitation of an engine is determined by the 
temperature at which any engine component would fail. Because of 
the poor response characteristics inherent in a system of tempera- 
ture control dependent on the sensing of metal temperatures, such a 
method is considered impracticable. The sensing and the control of 
combustion-gas temperatures appear to be more advantageous than an 
attempt to control directly the temperature of the critical engine 
component . 

Limiting control of gas temperatures presents a difficult 
problem in that direct means of sensing combustion-gas temperatures 
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of the magnitude encountered in gas-turbine engines are as yet 
unsatisfactory. The usual means of sensing temperature, such as 
thermocouples and resistance wires, are subject to errors resulting 
from radiation, conduction, and oxidation. In addition to these 
errors, the life of such instruments is comparatively short at the 
high temperatures to which they would be subjected. The thermo- 
dynamic methods presented herein were conceived to circumvent such 
difficulties. 

Any means of sensing temperature that is selected for control 
purposes should be accurate, respond quickly, and provide a response 
that is easily incorporated into a control system. Measurement and 
interpretation of the thermodynamic changes that occur in the work- 
ing fluid during the combustion process offer possibilities for 
meeting these requirements and obtaining temperature indications 
that can be used for control application. Temperature changes of 
the combustion gases are accompanied by pressure changes that are 
in effect instantaneous, are unaffected by errors of radiation and 
conduction, and can be utilized in a control system without modi- 
fication or amplification. 

Other investigations have been made correlating temperature 
rise across the combustion chamber with pressure changes. These 
correlations, however, depend on combustion- chamber pressure losses 
requiring a calibration of the particular combustion chamber and 
the method appears impractical for control purposes. 

A theoretical equation based on the thermodynamic properties 
of the combustion gas that correlates gas- temperature ratio with 
pressures upstream and downstream of the combustion zone has been 
developed at the NACA Cleveland laboratory and is presented 
herein. Means of simplifying this equation for control purposes 
are discussed and hypothetical control systems are presented to 
illustrate application of the control equations. 


DERIVATION OF CONTROL EQUATIONS 
General Control Equation 

Determination of the final temperature of the gas from the 
thermodynamic relations involved is possible without direct 
sensing of the temperature in the hot gas stream by equating the 
mass flow of gas in a duct ahead of the point of addition of heat 
(station 1) to the mass flow downstream of the point of addition 
of heat (station 2). 
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The symbols used herein are defined in appendix A. 

The flow at station 2 equals the air flow at station 1 plus 
the fuel flow f , that is, 

a 2 = a i + f 
a 2 = ( 1 + ij) a l 

When the equivalent gas-flow equations are substituted for a^ 
and &2 > 

* 2*2^ 1 | 28P 2 4P 2 - (* + Hi) *1*1^1 (1) 

P . 

The density p may be replaced by its equivalent — and p by its 

, ^ ^ 

corresponding relation 0' as shown in appendix B: 



The unknown variable,- total temperature Tg, may now be obtained 
from the equation 



Equation (3) offers a means of determining total temperature Tg 
and may be simplified for control application by using the following 
assumptions: 

1. The area ratio Ag/A^ is a known constant and may be 
obtained by direct measurement ; or by use of equation (3) it may be 
found by calibrating the system without heating the gas. 

2. In a particular design, the operating range of temperatures 
and pressure ratios is comparatively small and the conversion 
factor 0 may be selected as constant . The factor P is a multi- 
plication factor by which the hydraulic equation is converted to the 
compressible-flow equation. The equation for the conversion factor 0 
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therefore varies with each combination of static pressure p or 
total pressure F and static temperature t or total tempera- 
ture T used in the hydraulic equation so that the campressible- 
f low equation is always obtained. The expression for determining 

the actual value of the conversion factor 0' when p = ^ is 

developed in appendix B for the case in which the density p is 
proportional to the static pressure p divided by the total tem- 
perature T« This density relation is desirable for this analysis 
because it incorporates the total rather than the static tempera- 
ture and also because the numerical value of the conversion fac- 
tor 0' for each pressure ratio and value of the ratio of specific 
heats 7, where 7 is a function of temperature, is more nearly 
equal to 1 than for any other relation of the density p. The 
deviation of the conversion factor -0' from 1 is greatest at 
critical flow (fig. 1) and increases with Increasing temperature 
because of the change in the ratio of specific heats 7. At the 
point of critical flow and a value of 7 of 1.3, which corresponds 
to an air temperature of about 3000° R, the conversion factor 0' 
is approximately 0.945. The curve for 7 = 1.4, which represents 
an air temperature of about 500° E, is also shown. 

3. The area multiplier for thermal expansion Eq has a 
negligible change because relatively little change in temperature 
occurs at station 1 over the range of engine operating conditions. 


4. The area multiplier for thermal expansion Eg Increases 
to a value slightly above 1 with an increase in the fuel-air ratio 

E, 

remains a constant. For example, with 


f/a^ and the ratio 


1 + J- 


an assumed burner efficiency of 50 percent, the value of 


1+X 


a l 


decreases from 0.995 to 0.990 with a temperature differential 
increasing from 400° to 1600° F between stations 1 and 2. 


If the factors 
bined into a single 


t'zft'v ®1> ^ ppV 

a i 

constant K, equation (3) becomes 


are com- 


*2 


1 


(4) 
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This equation may "be adapted for control application. The average 
values of velocity pressures APi and AP2 may he obtained by 
using single pitot-static tubes at each station in a gas stream of 
known velocity distribution. 

Equation (4) may be further simplified for the case in which 
the ratio pg/p^ may be considered constant. Although an investi- 
gation should be made into the feasibility of making this assumption 
for any particular application, it will serve to simplify equa- 
tion (4) for illustrating a unique hydraulic control system subse- 
quently explained. For this condition, combination of the con- 
stants and Pg/p^ into a single constant K* and substitu- 

tion into equation (4) gives 

APg 

Tp = K' — - Tn (5) 

* APt 1 


Control Equation for Condition of Critical Flow 

Another simplification of temperature equation (3) may be made 
when critical flow exists at station 2 after the gas is heated. 
Equation (3) may be rewritten as - • 


Eo Ao i 



P2 P 2 ~ P 2 
Pl APj_ L 1 


( 6 ) 


where Pg - p 2 replaces its equivalent APg, and 0'g is ractored 
out of the bracketed term. 

For the condition of critical flow, the pressure ratio Pg/Pg 
is 




( 7 ) 


This equation may be substituted in equation (B8) (appendix B) for 
its equivalent r and an expression for (0'g)^ in terms of 7 is 
obtained 
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7-1 


■( M 1 


( 8 ) 


The terms in the factor 


E, 


2 *2 


.‘I* 1 


in equation (6) 


are assumed constant in the manner previously discussed for deter- 
mining the constant K in equation (4) and are designated by the 
constant B. Substituting B and the equivalent of pg from 
equation (7) and (0' 2 ) 2 from equation (8) causes equation (6) to 
become 

■" _1_ 

7-1 |_ / 2 \ 7_1 _ _ / 2 \7-l 


T 2 = B 2 


Mr) 


'-(vhfl 


(tIi 


M 


AP, 


This equation simplifies to 


T 2 = B2 |(f 


H) 


2+1 

7-1 


p 2 

r 2 


Pi AP, A 1 


(9) 


2+1 


The factor ^ ( y + l )^~^ var ies about 5 percent with a change 

in temperature from 500° to 3000° E. In the case of a temperature- 
limiting control for vhlch the value of total temperature T 2 is 

predetermined, the value of this factor may also be predetermined. 

A constant B' may then be used to replace the term 


B 2 2L (-*-) 
2 \ 7 + 1/ 

for a particular value of T 2 and equation (9) may be expressed as 


y+1 

7-1 
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B' 



T 

AP 


1 

1 


( 10 ) 


In the operation of current gas-turbine engines in which the 
ratio P 2 /Pi remains substantially constant, equation (10) may be 

written as 


B" 


Pi 2 Ti 

Pi AP X 


( 11 ) 


where B M is a constant. This control equation is unusual because 
sensing elements are required only before the combustion zone. 

The application of a control operating in accordance with 
equation (10) or (ll) can be illustrated by reference to the e 
state operating curve of a typical turbojet engine (fig. 2). The 
control operating line is shown superimposed on this engine- 
operating curve. In this engine, critical flow exists in the tur- 
bine nozzle box for engine speeds above 8000 rpm. For speeds less 
than this value, the nozzle-box Mach number M is less than 1. 

The maximum operating temperature shown, 2000° E, is the maximum 
permissible gas temperature specified for this particular engine. 
Below an engine speed of 8000 rpm, the temperature allowed by the 
control decreases with decreasing engine speeds; the values shown 
in figure 2 were determined from equation (10) and from typical 
engine data. The difference between the steady-state temperature 
at any engine speed and altitude and the corresponding temperature 
allowed by the control is the temperature differential available 
for acceleration. The maximum permissible gas temperature is 
allowed by the control only for those conditions of engine operation 
for which critical flow exists. Most normal operation occurs at 
engine speeds above 8000 rpm where the control allows operation at 
the maximum permissible temperature. From the shape of the 
operating curve, the temperatures required for steady-state opera- 
tion at low engine speeds apparently would not be allowed by the 
control. Making the control inoperative below some speed that 
must be determined from engine data therefore becomes necessary. 


APPLICATION OF CONTROL EQUATIONS 

A system designed to control total temperature T 2 by action 
based on the values of the other variables of equation (5) is 
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schematically shown in figure 3. Equation (5) was selected for 
illustration because a simpler presentation can be made than with 
the more basic and more generally applicable equation (4). This 
control system is unusual in that the required multiplication and 
division are accomplished by hydraulic means. Temperature control 
is maintained by the automatic regulation of resistance valve a 
in the fuel line, which regulates the fuel flow to the nozzles to 
limit the gas temperature to the predetermined maximum value . 

Valves b, c, and d can be located in any convenient fluid-flow 
line through which continuous flow is maintained, such as a fuel- 
bypass line or a lubricating-oil line. The following relation is 
obtained by equating the mass rate of flow through valves c and d: 


C c A c Vdp^ = 


or 




( 12 ) 


The temperature at station 1 So controls the valve area A^ that 
the following relation is maintained: 

= GTq 

where G is a constant. 

The velocity pressure AP^ is applied to diaphragm-operated 
valve b, which maintains the pressure drop across valve c Ap c 
equal to AP 1 . The velocity pressure AP 2 is applied across 
diaphragm-operated valve a, - which controls the fuel flow to the 
nozzles and thereby controls the magnitude of AP 2 . Valve a is 
in balance when the pressure drop across valve d Ap d equals AP 2 . 
In equation (12) therefore substitutions may be made for A^ 2 * 

Ap d , and Ap c : 



When equation (13) Is compared with equation (5), total tem- 

Ki /C „\ 2 * 


perature T 2 is found to equal 1 j A^. Area Aq may 
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therefore he set to give a predetermined value of - total tempera- 
ture Tg. If the actual value of total -temperature- T 2 is less 
than this value> the value of -velocity pressure AP 2 is less than 
static pressure difference Ap^ and valve a tends 'to open, 
allowing more fuel to enter the nozzles. If the actual value of 
total temperature T 2 is greater than this predetermined value, 
the opposite action occurs. If there are no other controls in the 
fuel-supply line to the nozzles, the control tends to maintain the 
total temperature T 2 as determined hy area Ac. If, however, 
controls exist in the system that restrict fuel flow in accordance 
with the d eman ds of requirements other than maintenance of a maxi- 
mum gas temperature, this control will act only to limit the maxi- 
mum value of total temperature Tg to the value determined hy the 
setting of area A c . 

A modification of this control is shown in figure 4. The 
operation of this second control is essentially the same au that of 
the first (fig. 3) except that a motion that is linear with changes 
in gas temperature is obtained hy the addition of a servo system. 
This temperature-dependent motion may he utilized in a control 
system that uses temperature as a control parameter in conjunction 
with other engine parameters for power-plant operation. Equa- 
tions (ll) and (12) are still valid, the only difference between 
the systems being in the method of control.. The static-pressure 
drop across valve c Ap c is controlled hy a diaphragm- operated 
servovalve and a servopiston to maintain static-pressure drop Ap c 
equal to the pressure difference AP 1 . Valve h is a diaphragm- 
operated throttling valve, which acts to maintain Ap^ equal to 
APg. The manual lever controls a servovalve and a servopiston- 
regulating valve a to maintain a set total temperature T 2 . 
Because area Ac is a measure of total temperature T 2 , the posi- 
tion of the stan of valve c is a measure of temperature as indi- 
cated hy the temperature indicator. The manual lever may he so 
set that the servovalve and the piston positioning valve a deter- 
mine the equilibrium position of the stem of valve c, which, in 
effect, determines the total temperature T 2 . 

Valves b, c, and d (fig. 4) in the bypass line are con- 
tinuously functioning and the position of valve c is a continuous 
indication of temperature. The temperature Indicator therefore 
operates when the control is used either as a temperature- limiting 
control or as a regulating control. This control system differs 
from the system shown in figure 3 in that the valve area A c 
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(fig. 3) is manually set by the pilot or may be preset to maintain 
a particular limiting temperature. This area setting indicates 
temperature only when the control is used as a regulating control; 
otherwise, the area setting is that of the limiting temperature. 

A control operating according to the control equations based 
on critical flow at one station is schematically shown in figure 5. 
Either equation (10) or (ll) maybe used as the temperature- control 
equation, depending on whether total pressure P 2 or P]_ is to 
be used to actuate the control. The schematic diagram in figure 5 
is based on equation (ll). This control is similar to those 
diagrammed in figures 3 and 4 in that valve a controls the fuel 
flow to the engine, and that valves b, c, and d are located 
in a convenient fuel-flow line . The density term p 1 /T 1 in equa- 
tion (ll) is so applied to a density-compensating system that 
valve b controls the' pressure drop Ap c proportionally to 
T 1 /p 1 « The area A c is made proportional to the total pressure P 1 
and the area A^ is made proportional to the constant B 2 of equa- 
tion (ll), which depends on the value of total temperature T 2 
selected. The static-pressure difference Ap^ then so controls 
valve a that AP^ equals Ap^,. When these substitutions are 
made, control equation (12) is seen to operate according to temper- 
ature equation (ll). 


CONCLUSIONS 

The analysis indicates that deteimination of combustion-gas 
temperatures from the thermodynamic relations involving gas temper- 
ature before combustion and gas pressures before and after combus- 
tion appears practical for high-temperature application such as for 
gas-turbine temperature control. For critical flow in the turbine 
nozzles, the analysis indicates that combustion-gas temperature can 
be determined and controlled from measurements taken only upstream 
of the combustion zone. 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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APPENDIX A 
SYMBOLS 

The following symbols were used in the report: 
area, sq ft 

air or gas flow, lb /sec 
coefficient of discharge 

area multiplier for thermal expansion of duct 
fuel flow, lb/sec 

acceleration due to gravity, ft /sec- 
Mach number 

total pressure, lb/sq ft absolute 
static pressure, lb/sq ft absolute 
total pressure minus static pressure, P-p 
gas constant 

static-pressure to total-pressure ratio, p/P 
total temperature, °E 
static temperature, °E 
velocity, ft/sec 

ratio of specific heats at constant pressure and constant 
volume 

arbitrarily selected density, lb/cu ft 
free-stream density, lb/cu ft 
stagnation density, lb/cu ft 
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0 conversion factor of hydraulic equation to compressible 

flow equation 

0* particular value of 0 when p = JL 

ET 

B,G,K constants 
Subscripts: 

a,b,c,d valves a, b, c and d in control systems 

1 station 1 (before heating) 
station 2 (after heating) 


2 
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APPENDIX B 


COMPRESSIBLE-FLOW EQUATIONS 

The hydraulic equation for incompressible flow may be multi- 
plied by an appropriate conversion factor 0 to obtain the exact 
equation for compressible flow. This expression for the conversion 
factor 0 may be derived from the compressible-flow equation by 
factoring out the hydraulic equation so that the remaining factor 
is the expression for the conversion factor 0. 

Bernoulli's theorem for compressible flow may be written as 

1 

-.2 

£) <B1> 

-U/J 

The weight-flow rate is 

a = Ap 0 V (B2) 



Substituting equation (Bl) in equation (B2) and replacing stagna- 

1 


tion density with the equivalent adiabatic relation 



gives 



The free-stream density Pq may be replaced by its equivalent 
and the equation simplified: 
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a * A 


**(t*t) 


Et 


P H) 


p 

H) 1 


a = A 




£ 

Et 


(?) 


1 


(B3) 


This egression for compressible flow may be written as 

a = A0 ^2gP (P - p) 


(B4) 


where 



2=1 ] 

i 1 ( 7 \ P 2 

/P\ 7 , 

p It - p) (rtrjRt 

\\3) ' 


(B5) 


This expression of 0 involves the density p, which also appears 
in equation (B4) and may be arbitrarily selected as a ratio 
involving total pressure P or static pressure p divided by 
total temperature T or static temperature t. For the case in 
which the density p is selected as p/RT, the conversion factor 0 
is designated as 0' and equation (B5) is simplified as follows: 


U >') 2 - 


Tp7 


, 2=1 

1(f)’ 


- 1 


■t (p - p)(r - i) 

The adiabatic relation of the temperatures is 


(B6) 


1 m (l\ 

t \P/ 


ill 

7 


(B7) 


This relation may be substituted into equation (B6) to obtain 
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If the pressure ratio p/P is set equal to r this equation 
becomes 




A plot of the conversion factor 0' against the pressure 
ratio r is presented in figure 1 , which shows the error that may 
be expected from neglecting 0 in the hydraulic equation (B4) 
where the density p is p/ET. The greatest deviation of the con- 
version factor 0 ' from 1 occurs at the critical pressure ratio, 
at which 0' is approximately 0.945 for the ratio of specific 
heats 7 equal to 1.3 for air at a temperature of 3000° B. 




total pressure 
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Pitot -static tube 



From fuel pump 



Figure 3. - Schematic diagram of basic temperature control 
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Pitot-static tube 



Figure 4. - Schematic diagram of temperature control with automatic temperature indication 
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